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We show a zero-gap semiconducting (ZGS) state in the quasi-two-dimensional organic con- 
ductor a-(BEDT-TTF)2l3 salt, which emerges under uniaxial pressure along the o-axis (the 
stacking axis of the BEDT-TTF molecule). The ZGS state is the state in which a Dirac cone 
with the band spectrum of a linear dispersion exists around the Fermi point connecting an unoc- 
cupied (electron) band with an occupied (hole) band. The spectrum exhibits a large anisotropy 
in velocity, which depends on the direction from the Fermi point. By varying the magnitude of 
several transfer energies of a tight-binding model with four sites per unit cell, it is shown that 
the ZGS state exists in a wide pressure range, and is attributable to the large anisotropy of the 
transfer energies along the stacking axis. 

KEYWORDS: organic conductor, BEDT-TTF, quarter-filling, tight-binding model, zero-gap semiconduc- 
tor, Dirac cone, contact point, a-(BEDT-TTF)2l3, 
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1. Introduction 

Quasi-two-dimensional organic conductors, BEDT- 
TTF (bis(ethylene-dithio)tetrathiafulvalene) salts, ex- 
hibit several electronic states, such as Mott insulator, 
charge ordering and superconductivity, 1 ' 2 under the vari- 
ation of temperature, hydrostatic pressure or uniaxial 
strain. The electronic states of these BEDT-TTF (ET) 
salts originate from 7r-electrons with the several trans- 
fer energies and the strong Coulomb interaction acting 
on-site and between nearest-neighbor sites. Among these 
salts, unusual phenomena have been observed in the a- 
(ET) 2 I 3 (a phase of BEDT-TTF tri-iodide) salt, 3 which 
consists of four molecules per unit cell with a 3/4- filled 
band (i.e., quarter-filled band of the hole). Under am- 
bient pressure (A. P.), the cv-(ET)2l3 salt exhibits the 
metal-insulator transition at Tmi = 135 K, which is fol- 
lowed by charge ordering along the 6-direction per- 
pendicular to the stacking axis (a-axis). When a uni- 
axial strain at 2 kbar is applied along the a-direction, 
the salt undergoes the superconducting transition at 7 
K. 10,11 Under high pressures, the a-(ET)2l3 salt shows 
exotic properties, that are consistent with the narrow 
gap semiconducting (NGS) state. 5,6 At 20 kbar of the 
hydrostatic pressure, resistivity becomes metallic and is 
almost independent of temperature between 300 K and 
1.5 K, although carrier density, n, decreases from 10 21 
cm~ 3 (300 K) to 10 15 cm" 3 (1 K). It has been claimed 
that such a temperature dependence may originate from 
the NGS state with a gap E g ~ 1 meV. 6 There exists a 
difference in the effect of uniaxial pressure between the 
a-direction (P a ) and the ^-direction (-P&). 11 For a com- 
pression along the a-direction, the pressure dependence 
of Tmi is small and the NGS state is likely to be real- 
ized for P a > 5 kbar. For the uniaxial strain along the 
6-direction, the metal-insulator transition is suppressed 
noticeably, whereas neither superconducting transition 
nor NGS state occurs at low temperatures. 

The ordered state in a-(ET)2l3 has been investigated 
theoretically; a charge ordering state with the horizontal 



pattern was obtained using the Hartree-Fock approxima- 
tion 12-15 and the superconducting state in the presence 
of the charge ordering was clarified in terms of the ran- 
dom phase approximation. 16, 17 It is noted that such a 
superconducting state occurs just before the onset of the 
NGS state under the a-axis pressure. However, the ori- 
gin of the NGS state is not yet clear, although only the 
uniaxial pressure along the a-axis leads to the NGS state 
in the tv-(ET) 2 i 3 salt. 

In the present paper, we examine the pressure depen- 
dence of the band structure of the a-(ET)2i3 salt within 
the tight-binding model in order to understand the NGS 
state, which emerges and remains in a wide range of the 
uniaxial pressure along the a-axis. In §2, we study the 
pressure dependence by assuming an extrapolation for- 
mula for the pressure dependence of transfer energies of 
the a-(ET)2l3 salt. We employ the data obtained from 
the X-ray diffraction experiment under pressures along 
both the a-axis and the 6-axis. 18-20 In §3, it is shown 
that the present organic conductor under high pressures 
exhibits a zero-gap semiconducting (ZGS) state with 
a Dirac-cone-like linear dispersion instead of the NGS 
state. The origin of the ZGS state under pressure is ana- 
lyzed by varying parameters of transfer energy in §4. We 
also investigate the ZGS state in §5 by using a reduced 
model with two sites in the unit cell. Section 6 is devoted 
to the summary and discussion. 

2. Formulation 

In order to describe conduction electrons for the a- 
(ET)2l3 salt, we start with a two dimensional lattice 
model, which consists of four ET molecules in the unit 
cell as shown in Fig. 1. The Hamiltonian is given by 



(1) 



where ti a: jp is the transfer energy from the (j,(3) site to 
the nearest-neighbor (n.n.) (i,a) site. In eq. (1), i and 
j (= 1, • • • , Nl) denote the unit cells of the square lat- 
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tice. The indices a and (3 denote the four molecules in 
the unit cell, where 1, 2, 3 and 4 correspond to A, A', B 
and C in ref. 18, respectively. q q is the annihilation op- 
erator of the conduction electron at the (i, a) site. Using 




Fig. 1. Conducting plane of ct-(ET)2l3, where there are four ET 
molecules in the unit cell and the bond represents eight transfer 
energies (t p i, ■ • ■ ,t p i,t c i, ■ ■ ■ ,t c 4). The z(j/)-axis corresponds to 
the 6(a)-axis in ref. 18. Thick arrows denote the dimerization in 
a model, which is considered in §5. 

the Fourier transformation, c ia — -/Vl -1 / 2 e lfc ' ri Cfc Q , 
cq. (1) is rewritten as 

H = £ a p{k)J ka c kp , (2) 

fe,a,/3 

where 

and tia-.jp represents eight transfer energies given by 
t p i,--- ,t P 4,t cl , ■■■ ,t ci as shown in Fig. 1. In cq. (2), 
s a p{k) is expressed as 



£aa (&) 


= (a= l,--- ,4) , 


e a p(k) 




£l2(fe) 


= td + t c2 e- lk y , 


£i3(fc) 




s 14 (k) 




£23(fc) 




£24(fc) 




e 34 (k) 


= t c3 + t ci e~ ik * , 



where k x is replaced by k x + 7r, and the suffix x{y) cor- 
responds to the 6(a)-axis. 18 The chemical potential de- 
termined by the 3/4 filling is set to zero in order to take 
the Fermi energy as the origin of the band spectrum. 

For the pressure dependence of transfer energies 
t ia :j0 = t A ,t' A (A = cl,- •• ,c4, pi,-- ■ ,pA) of a-(ET) 2 I 3 , 
we employ an extrapolation formula given by 16 

t A {P a ) = t A (l + K A P a ) , (5) 
t' A (P b ) = t A (l + K' A P b ) , (6) 



where t A (P a ) (t' A (Pb)) represents the transfer energy un- 
der the pressure P a (Pb) along the a-axis (6-axis). Here- 
after, the units of pressure and transfer energy are taken 
as kbar and eV, respectively. By using the data of a- 
(ET) 2 l3 salt for ambient pressure (P a = Pb = 0) 18 and 
those for P a = 2 and Pb = 3 at room temperature, 19,20 
the above coefficients are estimated as t p \ = 0.140, 
t p2 = 0.123, t p3 = -0.025, t pi = -0.062, t cl = 0.048, 
t C 2 = -0.020, t c3 = -0.028, t ci = -0.028 at ambient 
pressure and K p i = 0.011, K p2 = 0.0, K p3 = 0.0, K p 4 = 
0.032, K cl = 0.167, K c2 = -0.025, K c3 = 0.089 and 
K c4 = 0.089 {K' pl = 0.024, K' p2 = 0.031, K p3 = 0.053, 
K' pi = 0.022, K' cl = 0.042, K' c2 = 0.133, K' c3 = 0.167 
and K' c4 = 0.167) for pressures along the a-axis (6-axis). 
As seen from K A , the large K c i suggests that the pres- 
sure dependence along the a-axis is mainly determined 
by t c i. 

3. ZGS State 

In this section, we focus on the electronic states under 
the pressure along the a-axis, i.e., P a . Since there are 
four sites in the unit cell, there arc four energy bands 
which are obtained by diagonalizing cq. (2). The first 
band with the highest energy is almost empty and the 
second one is almost filled due to 3/4-filling. Whereas 
there is a small electron pocket close to the X(±7r, 0) 
point and a hole pocket close to the Y(0, ±7r) point at 
ambient pressure, 6 these pockets are reduced at P a = 
2i9,20 an( j are f ur ther reduced to a point in a wide 
range of higher pressures. 16, 17 Actually from eq. (5), it 
is found that such an unusual state is stable for P a > 3. 




Fig. 2. Band dispersion of a-(ET)2l3 at P a = 4 kbar in first 
Brillouin zone (left panel) and enlarged one around contact point 
at fe° = (0.602-7T, — 0.353-7r) (right panel), where Fermi energy is 
taken as origin. 



Figure 2 shows the energy band spectrum of a-(ET) 2 l3 
at P a = 4 (left panel) in the plane of k x and k y , which 
is scaled by the inverse of the lattice constant. The first 
band is unoccupied, whereas the second band is occu- 
pied. It should be noted that the first and second bands 
touch at the points ±fe° as shown with the circle. This 
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contact point becomes the Fermi point due to 3/4-filling. 
Actually, the matrix elements of eq. (4) , corresponding to 
fe° = (0.602vr, -0.353tt), are s 12 = 0.072 - 0.016i, e 13 = 
0.124 + 0.066i, £14 = 0.115 + 0.024i, e 23 = -0.135 - 
0.040i, e 2 4 = 0.014 - 0.117i and £ 34 = -0.055 - 0.034z. 
The four band energies are given as 0, 0, -0.210, -0.477, in 
which the Fermi energy is given by ep = 0.172. The ex- 
istence of such a point implies that the hole and electron 
bands are degenerate at the point. Since the point has 
no relation to the symmetry, such a degeneracy occurs 
accidentally as first pointed out by Herring. 21 

It should be noted in the enlarged figure (right panel 
of Fig. 2) that a linear dispersion is present around the 
contact point k . Such a state always exists for P a > 3. 
The dispersion close to the point can be expressed as 



E ± (Sk) = v ± (Sk)\k-k° 



(7) 



where 5k = k — k , and v + (Sk) (v~(8k)) is the velocity 
for the electron (hole). We call such a state the Dirac 
cone, according to the literature. 22-25 In Fig. 3, the ve- 
locity v (Sk) is shown as a function of 9, which is the 
angle between the vector Sk and the fc^-axis. A noticeable 
angular variation in velocity is seen with the relation: 

v+(6k) ^ -v-(5k) , (8) 

v+(Sk) = -v-(-Sk) , (9) 

within the numerical accuracy of the present calculation. 
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Fig. 3. Angular dependence of velocity around contact point, i.e., 
Dirac cone at P a =10 (solid line) and 5 (dashed line), where 9 
(horizontal axis) denotes angle between Sk(= k — k ) and the 
fee-axis. The inset shows (Sk) at P a = 10, where Sk > and 
< correspond to 8 = and n, respectively. 



For both low and high pressures, the schematic band 
dispersion is shown in Fig. 4, where the crossing in the 
band spectrum denotes the contact point which is located 
at the bottom of the electron band. At low pressures 




Fig. 4. Schematic band dispersion £ of metallic state at low pres- 
sure (a) and that of ZGS state at high pressure (b), where wave 
number k (horizontal axis) is chosen so as to include contact 
point. 



(Fig. 4(a)) with P a < 3, one finds a metallic state where 
two contact points at fc° and — k° exist below the Fermi 
energy ep. However, at high pressures (Fig. 4(b)) with 
P a > 3, the contact point is located just at the Fermi 
energy, i.e., the realization of the ZGS, which is retained 
due to the 3/4-filling for a wide pressure range. The ZGS 
state shows up in the density of states D(ui) close to the 
Fermi energy as 16 



D(u>) oc \lo\ 



(10) 



where u> = corresponds to the Fermi energy. Since D(to) 
vanishes at u = 0, we call the state with eq. (10) the ZGS 
state instead of the NGS state. 6 ' 16 Although the contact 
point still remains in the metallic state at low pressure, 
the corresponding cusp in the density of states is reduced 
and becomes invisible for P a < 2.5. 

In Fig. 5, the pressure dependence of the point contact, 
fe°, is examined, where the Fermi surface is shown for 
ambient pressure and P a = 2. With increasing pressure, 
the area of the Fermi surface is reduced as seen from the 
electron (hole) pocket around the X (Y) point. 19 The 
symbols on the arrow P a represent the location of contact 
points at P a = 0, 2 and 3, respectively. The arrow P a , 
which denotes the trajectory of k° from P a = to P a = 
10, extends to the direction of the T point with increasing 
P a ■ The contact point is located below the Fermi surface 
for P a < 3 and just on the Fermi surface for 3 < P a < 
10. The arrow Pf,, which corresponds to that of k° for 
< Pb < 10, is discussed in §6. 

4. Stability of Contact Point 

In order to examine the stability of the contact point 
against the change in parameters, we calculate the en- 
ergy band spectra by varying magnitudes of the transfer 
energies from those of a-(ET) 2 l3. Since the difference 
between t c \ and t c2 is large for a-(ET) 2 l3, we employ a 
model, in which t c \ ^ t c2 and the remaining parameters 
are chosen as 



pi 



p2 



1 , t 



p3 



Ep4 



(11) 



Thus, we examine a simplified model with independent 
parameters, t c \, t c2 , t p3 and t c3 . 

Figure 6 shows phase diagram in the plane of t c \ and 
t c2 for fixed t p3 (= t p 4 = t c3 = td) = —0.25. The con- 
tact points exist in regions I and II, which correspond to 
Figs. 4(b) and 4(a), respectively. Region I exhibits the 
ZGS state, where the point becomes equal to the Fermi 
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Contact Point 
AP. (0.660jc,-0.485ji) 
P a =2 (0.63171,-0.4 1 571) 
P a =3 (0.617tu,-0.383tc) 
P b =3 (0.671 tc,-0.407tc) 



Fig. 5. Fermi surface (FS) of a-(ET)2l3 at ambient pressure and 
P a = 2, where main (inset) figure corresponds to quarter (full) 
part of Brillouin zone. The arrow P a denotes the pressure depen- 
dence of the contact point k° for < P a < 10, whereas the arrow 
P b corresponds to fc° for < P„ < 10. The square, triangle, and 
filled circle represent kP for P a = 0, 2, and 3 respectively, while 
the open circle denotes that for P b = 3. 
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Fig. 6. Phase diagram on plane of t c \ and t c 2 for t v z = t p 4 = 
tcZ = tcA = —0.25, where transfer energies t c \ and t c 2 are also 
normalized by t p \{= t p 2). The arrow denotes the the variation in 
(tcli *c2) estimated from eq. (5) for < P a < 10 and < P D < 5, 
respectively, where the cross corresponds to ambient pressure. 



point. Region II corresponds to the metallic state with 
the point located below the Fermi energy. In Region III, 
the contact point disappears and the conventional metal- 
lic state is found. The boundary between regions II and 
III is given by 

t c i + t c2 = n, t c i + t c2 = r 2 , (12) 

where r\ = 1.53 and r 2 = —0.83. Such a straight line 
of the boundary is due to the choice of the parameter 
tp3 = tp4 = tc3 = t C 4- When parameters t c i and t c2 take 



the value on the boundary of r\ (r 2 ), the contact point 
is given by the T (X) point on the k x -k y plane in Fig. 5. 
Region I corresponds to a large t c i and a negative t c2 
(or a large t c2 and a negative t c i)- Note that the case of 
tci = t c2 (dotted line in Fig. 6) gives the metallic state 
with a special situation where the first band has the same 
energy as the second band on the line k y = ±7r in the 
Brillouin zone. The contact point disappears on this line. 
The van Hove singularity in the density of states exists 
near the Fermi energy when the parameters in region I 
move close to region II. When |£ P 3|(= |£ P 4|) and |i C 3|(= 
\t c i\) are decreased, the area of region II is reduced. 

Here, we comment on the pressure dependence of the 
contact point in Fig. 6. Here the arrow denotes the vari- 
ation of t c i and t c2 under the uniaxial pressures, P a and 
Pb. From the comparison of these two arrows, it is found 
that the increase in P a leads to the ZGS state but that 
in Pb is away from the ZGS state since the increase in 
t c \ is large enough for P a but is small for 

5. Reduced Model for ZGS State 

The condition for the ZGS state with contact points 
depends on the transfer energies, although some of them 
are redundant. In §3, the ZGS state was obtained for the 
model shown in Fig. 1. It consists of four sites in the unit 
cell. However, the model is complicated to examine the 
ZGS state analytically. In order to understand clearly the 
ZGS state, we reduce further the number of parameters 
to obtain a model with two sites per unit cell. 26 There 
are other reduced models reproducing the energy band 
spectrum of (ET) 2 X, 13 ' 15 which will be discussed in the 
next section. 

First, we ignore four transfer energies, t p z,t pi , t c3 and 
t C 4, in Fig. 1, since the ZGS state with contact points 
exists even for t p z — t p 4 = £ C 3 = t c ^ = 0. However, such 
a choice is not realistic for a-(ET) 2 l3 at P a = due to 
the absence of the metallic state as shown later. Next, for 
the simplicity of the calculation, we reduce the number of 
molecules by considering the dimerization for the bond 
t p2 between 1 and 4 sites (and also that for t p \ between 2 
and 4 sites), which is shown by the thick arrow in Fig. 1. 
With such a procedure, Fig. 1 can be replaced by Fig. 7, 
where we obtain an anisotropic square lattice consisting 
of two sites per unit cell. Thus, the reduced model is 
described by the Hamiltonian 

T 

T* 



H = 



t c ie 



t c2 e l 



(13) 



T — t p i + t p2 e 

with a half-filled band. Hereafter, we set t p i = t p2 (= 1) 
for simplicity. Two bands are degenerate at fc° when the 
off-diagonal component of eq. (13) vanishes. The wave 
vector, at which the degeneracy occurs, can be obtained 
analytically as 



k° = ±2 tan" 1 



fc° = ±2 tan" 



'r2 



2 + td + t c2 , 

^{2-tci-t c2 ){2 + t c i+tc 2 ) S 

t c2 — t c l 



(14) 
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Equation (13) gives a symmetric linear dispersion of the 
Dirac cone with v + (Sk) = —v~(Sk). It is slightly differ- 
ent from eq. (7). 




Fig. 7. Reduced model with anisotropic square lattice, which is 
obtained from Fig. 1. Here two sites are present in unit cell. 
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Fig. 8. Phase diagram on t c i-t c 2 plane for the reduced model 
given by Fig. 7. Region I is the same as that in Fig. 6, whereas 
region IV corresponds to the insulating state. The dotted line 
shows t c i = t c 2 , where the upper and lower bands are degenerate 
on the zone boundary. The inset shows the trajectory of two 
contact points, when t c \ and t c 2 are varied as the arrow in the 
main frame. 



Figure 8 shows the t cl -t c2 phase diagram obtained 
from the model of Fig. 7. The ZGS state is realized 
in region I and the insulating state is obtained in re- 
gion IV, respectively. The inset shows the locations of 
contact points in the k x -k y plane, which are obtained 
from cq. (14) when parameters t cl -t c2 are varied as 



(C t = -1,-2,-3) 

tc2=t cl +C t . (15) 

The corresponding parameters are shown by the arrow 
with 1, 2 and 3 in the main frame. The boundary between 
I and IV is given by 

t c2 = -t cl ±2 , (16) 

respectively. The properties of the contact point are as 
follows. The contact point is located at X for t c2 = —t c \ — 
2, and at T for t c2 = —t c i + 2, whereas the contact point 
is present near the zone boundary of the k x -k y plane for 
t c i ~ t c2 . When t c \ = t c2 holds, two bands are degenerate 
on the zone boundary and the contact point disappears. 
These facts indicate that the contact point comes from 
the symmetry breaking of the bond due to t c \ ^ t c2 . 

From Fig. 8, it is found that the ZGS state is mainly 
obtained for t c it c2 < 0, whereas the transfer energies 
{t P 3, tp4, t C 3, tci) are needed for the metallic state. 

6. Summary and Discussion 

In the present study, we examined the ZGS state 
by changing transfer energies within the tight-binding 
model. For a-(ET) 2 l3, the hole and electron bands are 
degenerate at two contact points k° and — k° in the first 
Brillouin zone. The dispersion of the electron and hole 
bands around the contact point resembles that of the 
massless Dirac Fermion, which is called the Dirac cone. 
In a wide pressure range (also parameters of transfer en- 
ergy range), there exists the ZGS state, in which the con- 
tact point coincides with the Fermi point. The location 
for the contact point moves in the k space continuously 
when transfer energies are varied. The stability of the 
contact point was examined using a simplified model. 
The contact point moves to the symmetry axis in the 
Brillouin zone when the parameters are varied toward 
those with higher symmetry. This was also verified using 
the reduced model. 

Let us comment on the state under the 6-axis pres- 
sure Pb, which is calculated from eq. (6). When p, is in- 
creased, the contact point still exists and moves as shown 
by the arrow P, (< 10) in Fig. 5. However, the ZGS state 
does not appear, since the contact point is always located 
below the Fermi energy. The absence of the ZGS state 
comes from the fact that the parameters t c \ and t c2 still 
remain in region II of Fig. 6. In other words, t c \ for p, 
does not increase much compared with that for P a . 

Let us examine if the ZGS state is present in other 
reduced models for a-(ET)2l3- Hotta 15 examined a dif- 
ferent model with two molecules in the unit cell by taking 
the dimer for the bond t p \ between 1 and 3 sites (and also 
that for t p2 between 2 and 4 sites) in Fig. 1, where the 
five transfer energies are kept. By substituting into these 
energies, the experimental data at ambient pressure, 20 
the contact point below the Fermi energy is obtained at 
A; = (0.6417T, — 0.69471"). This corresponds to the metal- 
lic state, although the location of k° differs from that 
of Fig. 2. For such a model, we also obtained the ZGS 
state under the pressure of P a . Comparing Hotta's model 
with the present one, both models take into account the 
dimerization for the bond with the large transfer energy, 
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but there is a difference in the choice of the bond as 
seen from Fig. 1. These two kinds of pairings, which are 
introduced for simplicity of the calculation, seem to be 
almost degenerate in their energy. In the present model 
of §5, the condition of t p3 = t p 4 = i c3 = t c4 = is 
further introduced in order to focus on the ZGS state. 
There is also another reduced model for a-(ET) 2 i3, 13 in 
which t c i = t c2 = t c3 = t c4 is taken. This model does not 
exhibit the ZGS state within the numerical calculation. 
This result may be ascribed to the fact that the ZGS 
state is expected mainly for t c it c2 < in Fig. 6. 

Those contact points obtained in the present study 
are related to the "accidental degeneracy" of the bands, 
which was pointed out by Herring. 21 In the case of "acci- 
dental degeneracy" , the degeneracy is not caused by the 
symmetry of the materials, but is caused by the intrinsic 
properties of the special band. As for the general case 
of the three-dimensional system, there appear contact 
lines in the Brillouin zone, which lead to a metal due to 
the existence of the Fermi surface. The ZGS state in a- 
(ET) 2 l3 exhibits the contact point, that always coincides 
with the Fermi energy, due to the quarter-filling in the 
highly two-dimensional system. The existence of such a 
degeneracy in a-(ET) 2 i3 has been supported recently by 
first-principle calculations. 27 ' 28 

Anomalous properties associated with the Dirac cone 
have been investigated for bismuth 22,23 and graphite 24, 29 
using the effective Hamiltonian for the Dirac cone. In 
those cases the degeneracy occurs on the symmetry axis 
of its Brillouin zone. The effect of Berry phase due to 
the Dirac cone has been observed using the phase anal- 
ysis of quantum oscillations. 25 Thus, the ZGS state of 
a-(ET) 2 I 3 is also expected to exhibit many anomalous 
properties. The new phenomenon known as anomalous 
transport 5 ' 6,30 may come from the Dirac cone and the 
unique properties of the ZGS state. 
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